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I .   INTRODUCTION 

The  electro-optic  effect  is  an  interaction  between  high-frequency 
electromagnetic  fields  (such  as  light)  and  low- frequency  electric  fields 
(such  as  static  fields  or  microwaves).   The  electro-optic  effect  is  an  example 
of  a  large  class  of  phenomena  (including  optical  second-harmonic  generation, 
Stoke 's  radiation,  and  the  stimulated  Raman  effect)  generally  called  "nonlinear 
optics."   '      Certain  materials,  such  as  the  liquid  nitrobenzene  and  the 


crystals  quartz  and  KH  PO,  (KDP) ,  exhibit  a  strong  (measurable)  electro-optic 
effect,         as  well  as  other  nonlinear  optical  effects. 

For  analysis,  the  electro-optic  effect  is  well  approximated  by 
expressing  the  optical  dielectric  properties  of  the  material  as  a  function  of 
the  local  strength  of  the  applied  static  field.   Thus,  for  example,  the  local 
phase  velocity  (or  index  of  refraction)  becomes  a  function  of  the  static 
field,  so  that  the  relative  phase  retardation  of  light  passing  through  the 
material  may  be  controlled  by  externally  applied  static  or  radio- frequency 
fields . 

In  application  the  electro-optic  effect  is  the  basis  for  a  variety 
of  devices  called  "electro-optic  light  modulators."  Perhaps  the  most  familiar 
of  these  devices  is  the  Kerr  cell  (or  "light  shutter"),  used  since  the  1920's 
for  high-speed  photography.   In  recent  years  the  electro-optic  effect  has 


been  used  in  more  sophisticated  devices  to  produce  modulation  of  phase,  polar 

of  light, '•3^' 
W,    [5],  [6] 


ization,  and  amplitude  of  light,    '     and  to  control  the  direction  of 


light  beam  propagation. 

Most  electro-optic  light  modulators  are  designed  for  time  modulation 
of  light,  similar  to  modulation  of  a  radio-frequency  carrier.   The  "electro- 
optic  light  valve"  (sometimes  called  an  "Ardenne  tube")  is  a  spatial  modulation 


.ce;  that  is,  the  light  valve  imposes  a  modulation  pattern  'across  the 
surface  area  of  a  wavefront,  analogous  to  the  intensity-modulated  wavefront 
of  light  which  has  passed  through  a  photographic  transparency.   The  light 
valve  is  potentially  capable  of  reproducing  a  television  image  in  the  form  of 
a  "real-time  photographic  transparency. "  This  potential  has  resulted  in 

investigations  in  the  use  of  the  electro-optic  light  valve  as  part  of  a 

[71 
television  projection  system    and  as  the  image  source  for  real-time  optical 

data  processing. 

The  configuration  of  a  typical  electro-optic  light  valve  is  shown 

in  Figure  1.   The  light  valve  consists  of  a  vacuum  chamber  containing  an 

electron  gun  and  a  large  flat  thin  electro-optic  crystal,  such  as  KDP.   One 

face  of  the  crystal  is  covered  by  a  transparent  conducting  coating  which  is 

electrically  grounded.   The  electron  gun  deposits  a  charge  pattern  on  the 

other  face  of  the  crystal.   The  charge  pattern  induces  a  static  electric 

field  inside  the  crystal,  so  that  light  passing  through  is  phase-modulated 

by  the  electro-optic  effect.   If  the  electron  beam  is  scanned  over  the  crystal 

face  in  a  raster  pattern  and  the  beam  current  is  controlled  by  a  television 

video  signal,  the  resulting  charge  pattern  is  a  representation  of  a  television 

picture.   If  the  crystal  is  thin  enough,  the  local  charge  density  controls  the 

strength  of  the  local  electric  field,  which  in  turn  controls  the  local  phase 

retardation  of  the  emerging  wavefront  of  light.   Thus  the  emerging  light  is 

spatially  modulated  into  a  pattern  resembling  the  original  television  image. 

. sible  (intensity  modulated)  image  may  be  obtained  if  the  light  is  polarized 

before  entering  the  crystal,  and  passes  through  an  analyzer  after  emerging 

from  the  crystal  (as  in  Figure  1). 
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Various  versions  of  the  electro-optic  light  valve  have  been  under 
development  for  about  ten  years.  The  physical  con- 

struction and  maintenance  of  these  devices  have  presented  many  difficulties, 
but  most  have  been  solved  by  developing  suitable  technologies .  A  more  serious 
problem  has  yet  to  be  resolved:  the  modulation  characteristics  have  been 

:locre  at  best,  so  that  -  in  particular  -  the  maximum  pattern  resolution 
currently  feasible  is  less  than  25%  of  standard  television  resolution. 

As  will  be  seen,  the  pattern  resolution  is  primarily  determined  by 

the  thickness  of  the  crystal  relative  to  the  size  of  its  faces.   This  problem 

[9] 
was  recognized  quite  early,    and  various  attempts  have  been  made  to  develop 

a  theoretical  analysis  of  the  light  valve  modulation  characteristics  in  order 
to  provide  a  quantitative  explanation  of  resolution.  We  now  briefly  consider 
two  of  these  analyses . 

The  simplest  (and  perhaps  earliest)  analysis    considers  a  disk 
of  diamater  D  located  a  distance  T  from  a  grounded  conducting  plane.   The 
disk  is  charged  to  a  potential  V,  and  in  the  plane  of  the  disk  the  surrounding 
potential  falls  to  approximately  V/2  at  a  diameter  D  +  O.M+T.   Assuming  that 
the  intensity  (power)  of  the  transmitted  light  is  proportional  to  the  potential 
in  the  plane  of  the  spot,  the  half -power  diameter  of  the  light  spot  is  equal 
to  the  V/2  diameter  above.   In  general  the  electron  beam  spot  size  (~  D)  is 
smaller  than  T,  so  the  half-power  diameter  is  approximately  T/2.   If  we  define 
the  minimum  resolvable  separation  of  two  light  spots  as  approximately  twice 
the  half-power  diameter,  then  resolution  spacing  is  roughly  T,  the  thickness 

the  crystal.  Standard  television  resolution  is  500  line-pairs,  correspon- 
ding to  500  spots  in  a  line,  so  that  this  resolution  would  require  that  the 
face  of  the  crystal  be  at  least  500  times  the  crystal  thickness.  Thus  a  1" 
diameter  crystal  must  be  polished  to  less  than  2  mils  in  thickness. 


[12]   [131 
The  second  analysis    '      is  considerably' more  detailed  and 

explicit  than  the  first,  but  follows  similar  arguments.   Again  we  have  a  disk 

at  potential  V  located  a  distance  T  from  the  grounded  conducting  plane,  but 

a  quantitative  expression  is  derived  for  the  potential  at  all  points  in  the 

plane  of  the  disk.   Further,  it  is  assumed  that  the  intensity  of  transmitted 

light  is  proportional  to  sin  [a   V(p)  +  0],  -where  V(p)  is  the  disk-plane 

potential  at  a  distance  p  from  the  spot  center  and  where  a  and  p  are  suitable 

2 
constants  (this  sin  dependence  is  found  in  typical  electro-optic  light 

[3] 

modulators    which  use  an  analyzer  to  convert  polarization  modulation  to 

intensity  modulation) .  The  resulting  expressions  for  the  intensity  distri- 
bution of  the  light  spot  are  analyzed  by  a  computer  program  to  obtain  reso- 
lution characteristics  for  a  set  of  charge  spot  diameters  and  crystal  thickness 

[13] 
Another  approach     is  used  to  obtain  the  normalized  Fourier  decomposition 

of  the  disk-surface  potential  V(x).   This  Fourier  decomposition  represents 

the  spatial  frequency  response  of  the  surface  potential  in  terms  of  the  charge 

pattern  and  can  be  used  to  obtain  an  equivalent  "pass-band"  for  the  device. 

This  (normalized)  Fourier  decomposition  V(N)/V(0)  is  given  by 

V(N)   1  -  e"2ltNT   r,    _x    *  nN  -2jtNT,-1 

vM  =   nKm         [(K  +  1)  +  (K  -  1)e    ] 

where  N  is  the  resolution  in  television  lines  (half -cycles)  per  millimeter, 
T  is  the  crystal  thickness  in  millimeters,  and  k  is  the  relative  dielectric 
constant  of  the  crystal.   This  spatial  response  analysis  gives  resolution 
characteristics  similar  to  those  of  the  light-spot  analysis.   The  analysis 
concludes  with  a  statement  predicting  that  for  a  crystal  of  thickness  T  = 
O.36  mm  and  a  resolution  of  500  television  lines,  the  crystal  face  must  be 
approximately  118  mm  square.   By  linear  interpolation  we  conclude  that  a  1" 
diameter  crystal  must  have  a  thickness  of  less  than  3  mils. 


The  quantitative  results  of  these  two  analyses  are  similar:   a  1" 
crystal  must  be  polished  to  a  thickness  of  2-3  mils.  Electro-optic  crystals 
of  useable  quality  (e.g.,  uniformity)  are  currently  available  in  1"  -  2" 
sizes,  but  crystals  of  this  size  have  not  yet  been  polished  to  a  thickness  of 
less  than  about  10  mils.   Thus  tolerable  resolution  from  an  electro-optic 
light  valve  is  in  essence  a  technological  problem,  at  least  until  crystals 
are  available  in  the  dimensions  indicated  by  the  preceding  two  analyses. 

In  considering  these  two  resolution  analyses,  one  finds  they  share 
two  assumptions:   First,  the  static  fields  are  determined  for  an  isotropic 
crystal;  and  second,  the  transmitted  light  intensity  is  considered  a  function 
of  the  surface  potential  only.   As  will  be  seen,  the  first  assumption  —  an 
isotropic  crystal  --  is  reasonable  since  even  though  the  electro-optic 
crystals  are  in  general  anisotropic,  the  static  field  equations  and  solutions 
have  the  same  functional  form  in  either  an  isotropic  or  an  anisotropic 
dielectric.   However,  the  resulting  quantitative  descriptions  are  somewhat 
different.  Tne   second  assumption  essentially  implies  that  each  element  of 
the  crystal  face  is  considered  as  a  local  intensity  modulator,  with  a  uniform 
static  field  in  the  direction  of  light  propagation  determined  by  the  local 
surface  potential  (as  in  a  parallel  plate  capacitor).   This  local  modulator 

is  analyzed  as  a  "standard"  longitudinal  (r^-)  linear  electro-optic  polari- 

r  3i 

zation  modulator.     That  is,  the  crystal  is  normally  uniaxial  but  becomes 

biaxial  in  the  presence  of  a  static  field  along  the  optic  (x  )  axis.   Light 
propagating  through  the  crystal  in  the  x  direction  splits  into  two  orthogonal 
components,  each  traveling  at  a  different  phase  velocity  controlled  by  the 
value  of  the  static  field.  Leaving  the  crystal  the  two  components  recombine 
with  a  new  polarization  determined  by  the  differential  phase  modulation  and 


thus  by  the  value  of  the  static  field.   As  a  first-order  approximation  for 
spatial  as  well  as  polarization  modulation  this  analysis  is  only  partially  valid, 
Spatial  modulation  requires  a  local  variation  in  the  direction  of  propagation, 
so  the  light  no  longer  travels  along  the  optic  axis.   Since  the  two  biaxial 
phase  velocities  are  a  function  of  direction  of  propagation,  we  might  expect 
that  pattern  direction  (in  addition  to  static-field  characteristics)  will 
determine  the  spatial  resolution  of  the  light  modulator.  As  will  be  seen  in 
the  first-order  analysis  that  follows,  the  resolution  characteristics  are  a 
function  of  direction,  at  least  in  theory. 

The  present  work  is  an  attempt  to  develop  a  more  rigorous 
theoretical  analysis  of  the  modulation  characteristics  of  the  electro-optic 
light  valve.   In  Section  II  the  complete  static  field  distribution  is  obtained 
for  an  anisotropic  crystal,  and  for  the  spatial  Fourier  decomposition  of  an 
arbitrary  charge  pattern.   This  results  in  the  anisotropic  equivalent  of  the 
resolution  characteristics  V(N)/V(0)  discussed  previously.   In  Section  III  the 
static  field  solution  is  combined  with  the  linear  electro-optic  dielectric 
tensor  and  Maxwell's  equations,  and  a  first-order  solution  of  the  resulting 
wave  equation  is  obtained  by  a.  perturbation  analysis.   This  first  order 
solution  is  in  the  form  of  a  far-field  radiation  pattern,  giving  a  direct 
expression  of  the  first-order  optical  modulation  characteristics  of  the 
electro-optic  light  valve . 

The  static  field  solution  and  the  first  order  optical  solution  are 
obtained  by  a  generalized  method  for  boundary  conditions  based  on  classical 
image  sources  (essentially  an  application  of  Green's  theorem).   These 
solutions  and  the  method  for  their  development  are,  to  the  best  of  the 
author ' s  knowledge ,  original . 


II.   ANISOTROPIC  STATIC  FIELD  SOLUTION 

The  static  field  equations  for  anisotropic  media  are  very  similar 
to  those  for  ordinary  isotropic  media,  since  here  the  material  is  assumed  to 
be  homogeneous.  However,  the  anisotropic  parameters  (the  dielectric  tensor 
c)  complicate  the  analysis  somewhat.   The  physical  configuration  of  the  light 
valve  is  such  that  the  general  static  field  equations  and  boundary  conditions 
can  be  simplified  considerably,  so  that  we  will  first  define  the  configuration 
to  be  analyzed. 

The  configuration  which  describes  the  static  field  conditions 
appropriate  for  the  light  valve  is  shown  in  Figure  2.   As  discussed  in  Section 
I,  the  static  field  is  established  by  a  surface  charge  density  located  on  one 
face  of  a  flat  thin  crystal.   This  crystal  has  been  oriented  and  cut  so  that 
the  surfaces  of  the  broad  flat  parallel  faces  are  normal  to  one  of  the 
principal  axes  of  the  crystal.   Let  the  coordinate  system  of  the  fields 
coincide  with  the  principal  axes  of  the  crystal  with,  say,  the  x  -  axis 
normal  to  the  crystal  faces.   Let  the  x  -  thickness  of  the  crystal  be  T  and 
choose  the  origin  and  direction  of  the  x  axis  so  that  the  grounded  con- 
ducting surface  is  on  the  right,  at  x  =0,  and  the  charged  surface  on  the 
left,  at  x0   =  -T.  Denote  the  surface  charge  distribution  on  the  left  by 
o(x, ,  x  )  (this  charge  layer  has  been  deposited  by  the  electron  beam).  Since 
the  size  of  the  charge  distributions  to  be  considered  is  small  compared  with 
that  of  the  face  of  the  crystal,  we  can  ignore  edge  effects  and  extend  the 
crystal  to  infinity  in  the  x, ,  x  directions.  Finally,  we  will  refer  to  the 
left  half-space,  x  <  -T,  as  Region  1  (free  space),  and  to  the  crystal  region, 

0,  as  Region  2. 


Region  1 
(Free  Space 
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^3 


Figure  2.   Crystal  Configuration  for  Static  Fields. 
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Maxwell's  equations  are  independent  of  the  properties  of  the 
material  so  that  as  usual  the  static  field  equations  are 

V  x  E  =0  (2.1) 

V  •  D  =  p  (2.2) 

By  (2.1)  we  may  express  the  electric  field  E  as  the  negative  gradient  of  a 
scalar  potential  $: 

E  =  -^J>  (2.3) 

For  a  linear  anisotropic  material  the  electric  material  equation  is  in 
general  a  second-order  tensor  relation  between  D  and  E.   Here  we  assume  the 
crystal  is  homogeneous  so  the  elements  of  the  relative  dielectric  tensor  e 
are  constants.   In  the  principal-axes  coordinate  system  e  is  diagonal,  and 
the  electric  material  equation  reduces  to 


D.  =  e  e.E.,  i  =  1,  2,  3  (2.4) 

l    o  1  i'      '   ' 


Using  (2.4),  the  divergence  equation  (2.2)  becomes 

I  >1E1  +  €2a2E2  +  €3a3E3  =  p/€o  (2*5) 

where  the  operator  d.  denotes  d/dx. .  In  terms  of  the  potential  $  the  above 
equation  (2-5)  is 


Ul\     +  e2^22  +  €3a32)  $  =  "p/€o  (2*6) 


11 


This  is  the  anisotropic  Poisson's  equation,  similar  to  the  isotropic  Poisson's 

equation  except  for  the  coefficients  e. .   In  fact,  if  we  use  the  coordinate 

transformation,   x.  =  u./FT,     then  (2.6)  in  the  u.  -  coordinate  system  is 
'111'  1 

an  isotropic  Poisson's  equation.  By  use  of  the  inverse  of  the  above  trans- 
formation, any  known  isotropic  solution  for  Poisson's  equation  can  be  con- 
verted into  the  equivalent  anisotropic  solution,  if  one  is  careful  to  observe 
that  this  coordinate  transformation  is  not  measure  -  preserving  (that  is, 
one  must  follow  the  general  tensor  rules  for  a  coordinate  transformation). 
Let  us  denote  the  potential  inside  the  crystal  (Region  2)  by  §  , 

and  that  in  free  space  (Region  l)  by  $  .   The  $  must  satisfy  (2.6)  with  the 

L         K 

e.  determined  by  the  particular  crystal  in  use.  However,  $  satisfies  the 
free-space  (isotropic)  Poisson's  equation,  which  is  (2.6)  .with  e  =  e  =  e  =  1, 

In  addition,  $  and  $  must  satisfy  the  boundary  conditions  at  the  inter- 

L       R 

faces  x_  =  0  and  x  =  -T.   The  general  anisotropic  boundary  conditions  are 

Ilk] 
rather  complicated,     but  in  this  case  the  surface  normals  coincide  with 

the  principal  axes  of  the  crystal.   Then  at  the  left-hand  interface,  con- 
tinuity of  the  normal  component  of  D  requires 


a30L  -  e333$R  =  c/e0  at  x3  =  -T  (2.7) 


where  a    is  the  surface  charge  distribution  described  above.   Continuity  of 
the  transverse  component  of  E  requires,  as  usual, 


°L  =  °R  at  X3  =  "T  (2*8) 


12 


Since  the  right-hand  interface  is  grounded  we  must  have 


3>R  =  0  at  x  =  0  (2.9) 


Finally,  since  this  grounded  plane  is  conducting,  the  transverse  component 
of  E  vanishes,  requiring 


\$R  =  ^2°R  =  °  at  x3  =  °  (2-10) 


These  boundary  conditions  (2.7)  through  (2.10),  together  with  the  appropriate 
forms  of  Poisson's  equation  (2.6),  are  necessary  and  sufficient  to  completely 
define  the  static  fields  in  terms  of  arbitrary  source  charge  distributions 
0   and  p.  Actually,  we  will  assume  that  p  vanishes  (since  the  resistivity  of 
the  crystal  is  very  large),  but  it  will  prove  convenient  to  retain  p  in  (2.6) 
until  later. 

We  will  obtain  the  static  field  solutions  by  using  Fourier  trans- 
forms. This  approach  is  advantageous  for  two  reasons:  first,  to  obtain  the 
optical  solutions  we  will  use  only  the  Fourier  transform  of  the  static  field 
solutions;  and  second,  in  this  particular  case  the  Fourier  transform  solution 
is  quite  simple.  We  assume  most  functions  of  interest  have  three-dimension 
Fourier  transforms,  so  that  some  scalar  function  f(x)  and  its  Fourier  trans- 
form f (k)  are  related  by 


f(k)  =  /*"  d3x  f (x)  e_tj-  '  -  (2.11) 

^ -00 

f(x)  -fCOJ^j3     f(k)  e*~   "  2  (2.12) 


13 


where 


k  •  x  =  k^  +  kgx2  +  k  x3 


By  differentiating  (2.12)  with  respect  to  x. ,  it  is  easy  to  show  that  the 
functions 


d.f(x),  Vf(x),  and  V^x) 


1  — 


have  as  their  respective  Fourier  transforms 


jk.f(k),  jkf(k),  and  -k2f(k) 


where 


,2   ,2   ,2   .2 
1  +  k2     3 


Thus  the  anisotropic  Poisson's  equation  (2.6),  which  is  in  x  -  space 
representation,  transforms  into  the  k  -  space  representation 


(elkl2  +  €2k22  +  e3k32)  0(k)  =  P&Veo  (2'13) 


while  the  k  -  space  isotropic  Poisson's  equation  is  simply 


k2  $(k)  =  p(k)/e  (2.U+) 


11+ 


The  Fourier  transforms  of  the  boundary  conditions  (two  -  dimension  trans- 
forms, since  x  is  fixed)  will  be  considered  later. 

Let  us  for  the  moment  consider  the  isotropic  Poisson's  equation 
in  k  -  space,  (2.l4).  For  a  unit  point  -  charge  located  at  x  =  a, 


p(x)  =  o(x  -  a)  =  6(x1  -  ax)o(x2  -  a2)o(x  -  a^) 
By  (2.11)  we  evidently  have 


/.  >.    -jk  •  a 
p(k)  =  e  °-      - 


so  that  $(k)  is  simply,  from  (2.l4), 


4>(k)  =  e"^  '  ^/eok2 


It  can  easily  be  shown  that  the  inverse  transform  of  $  above  is  given  by 

$(x)  =  1/kne   r 

where 

r  =  |x  -  a| 

This  $(x)  is  obviously  the  familiar  Coulomb  potential  of  a  unit  point  - 
charge.   However,  for  our  purposes  we  will  transform  $(k)  back  along  only 
one  direction,  say  x  : 
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(DC^,  k2>  x3)  =  J  2T    *(k)e    3*3      * 

J  2a€o  k/  +  k/  +  k32 


If  we  treat  this  integral  as  a  contour  integral,  there  are  two  poles  in  the 
complex  k  -  plane : 


k3  =  +  jTk^Tk^  =  +  jkr,  say  (2.15) 


This  is  an  ordinary  Fourier  integral  so  that  for  convergence  of  a  half  - 
plane  contour  there  are  two  cases  of  interest:   for  x   >  a   we  must  use 
an  upper  half  -  plane  contour  (lm(k  )  >  0),  and  for  x  <  a  we  must  use  a 

0  ~j  -J 

lower  half  -  plane  contour.   Then 


$(k  ,  k  ,  x  )  =  +  2njZ  (residues) 

-k  |x  -  a 
=  e-j(kiai  +  k2a2}  -$■ 


2k  e 
r  o 


This  result  is  not  remarkable  in  itself,  but  it  suggests  that  a  "mixed  - 
space"  representation  (that  is,  (k, ,  k  ,  x  )  -  space)  might  be  particularly 
useful  if  the  boundary  conditions  apply  to  x  only,  since  <3>(k  ,  k  ,  x  )  for 
a  point  charge  is  a  simple  exponential.  We  will  make  frequent  use  of  this 
mixed  -  space  representation  in  -what  follows. 

Next  we  consider  the  method  of  images.  Recall  the  elementary 
image  -  source  distribution  for  a  grounded  conducting  plane:   the  fields  on 
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one  side  of  the  plane  may  be  determined  by  replacing  the  conducting  plane 
and  opposite  half  -  space  with  a  negative  mirror  image  of  the  charge  dis- 
tribution in  the  half  -  space  of  interest.   In  our  case  the  conducting  plane 
is  defined  by  x  =  0.  Using  the  isotropic  Poisson's  equation  as  an  example, 
the  potential  $(x)  then  satisfies 

-V^Cx)  =  —  [pUl5  x  ,  x  )  -  p(x  ,  x  ,  -  x  )]  (2.16) 

o  ^ 

where  p(x)  is  assumed  to  vanish  for  x   >  0.   It  is  easily  shown  that  $(x) 
determined  by  (2.16)  completely  satisfies  the  required  boundary  conditions 
at  the  grounded  conducting  plane  ((2.9)  and  (2.10))  and  satisfies  the 
original  Poisson's  equation  using  only  p(x).   By  the  uniqueness  theorem  this 
is  the  solution  (in  the  half  -  space  x  <  0  in  this  example).   For  use  later, 
the  k  -  space  representation  of  (2.16)  is 


eok2<D(k)  =  p(kr  k2,  k3)  -  p(kr  kg,  -k3)  (2.17) 


since  from  (2.11) 


f(k.)  =  -f(-k.)  if  f(x.)  =  -f(-x.) 
x  i7        i'       i' 


The  method  of  images  may  be  generalized  to  analyze  a  variety  of 
physical  configurations.   As  an  example,  consider  an  object  -  say,  a  chess- 
man -  seen  through  a  thin  lens:   the  image  observed  is  a  magnified  copy  of 
the  original  object.   If  the  object  and  lens  are  replaced  by  an  exact  copy 
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of  the  object  magnified  to  the  proper  image  dimensions, -the  observed  image 
is  unchanged.   Alternately,  a  hologram  which  has  recorded  the  wavefront 
emerging  from  the  lens  will  also  reproduce  the  image  seen  by  the  observer. 
Both  the  hologram  and  the  magnified  copy  of  the  object  may  be  considered  as 
image  sources  in  the  sense  that  each  reproduces  the  observed  image  obtained 
from  the  original  configuration.  Note  that  in  these  two  cases  all  conditions 
in  the  observation  space  ("real"  space)  are  unchanged  even  though  conditions 
in  the  image  space  (all  space  outside  the  observation  space)  are  quite 
different.   This  example  of  generalized  image  sources  suggests  that  image 
space  conditions  on  the  image  source  are  ambiguous :  we  are  free  to  define 
image  space  conditions  in  any  self  -  consistent  manner  which  does  not 
disturb  the  original  image  in  observation  space.   This  statement  may  be 

ri^i     ri51 

rigorously  demonstrated  as  a  consequence  of  Green's  integral  theorems,    ' 

but  such  a  justification  is  beyond  the  scope  of  the  present  work.   Finally, 

it  is  worth  noting  that  the  image  source  concept  applies  not  only  to  true 

sources  located  outside  the  observation  space,  but  also  to  the  "reflected 

image"  caused  by  true  sources  inside  the  observation  space. 

Since  we  are  free  to  define  image  space  conditions,  let  us 

separate  the  static  field  problem  into  two  image  source  problems:   first,  a 

general  solution  for  the  potential  $  in  terms  of  an  image  source  p  located 

Li  h 

outside  Region  1,  where  §  and  p  satisfy  the  free  -  space  Poisson's  equation 
everywhere ;  and  second,  a  general  solution  for  $  in  terms  of  an  image  source 
p  located  outside  Region  2,  where  $  and  p^  everywhere  satisfy  the  aniso- 
tropic  Poisson's  equation.   The  final  solution  will  be  determined  by  requiring 
that  $  and  $  satisfy  the  boundary  conditions.   By  considering  two  separate 
full  spaces,  each  homogeneous  everywhere,  we  are  allowed  to  use  the  k  -  space 
Poisson's  equation  as  presented  above  (an  alternate  approach  would  be  to  use 
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the  mixed  -  space  Poisson's  equation,  which  in  this  case  would  be  a  piece  - 
wise  linear  ordinary  differential  equation). 

Accordingly,  from  (2.1U),  the  free  -  space  solution  is  determined 
by 

<DL(k)  =  pL(k)/eok2  (2.18) 

To  evaluate  the  boundary  conditions  we  must  use  the  mixed  -  space  represen- 
tation for  0T ,  obtained  from  the  ordinary  Fourier  integral 
h 

jk  x 
n    dk     p  (k)  e  *   J 

Vkl'  k2>  X3)  =  J4      k2  +  k2  &'& 

3     r 

2     2     2 
where  k  =  k0  +  k  "  as  in  (2.15).   We  are  interested  in  Region  1,  x„  <  -T, 
3    r  '  '  3 

and  by  definition  the  image  source  p  must  vanish  everywhere  in  Region  1, 

so  that  singularities  of  p  (k)  do  not  contribute  to  the  contour  integral  of 

Li 

(2.19)  evaluated  for  x  <  -T.   Then  the  only  singularity  of  interest  is  the 
lower  -  half  k  -  plane  pole  k  =  -jk  .  We  then  have,  from  (2.19), 

k  x_ 
r  3 

d>L(k1,  k2,  x  )  =  pL(k1,  k2)  1^^-   for  x  <  -T         (2.20) 

o  r 


where  PL(\,   ^2)  =  PL(k_L,  k2,  k^  =  -jk^) 


To  obtain  the  general  solution  for  <£>_  in  Region  2,  we  impose  two 

K 

conditions  on  the  image  source  pD.   First,  the  image  condition  requires  that 

K 

p  vanish  in  Region  2,  for  -T  <  x  <  0.   For  convenience  we  will  use  the 

jk  T 
form  p-.(x  +  T),  so  that  in  k  -  space  we  have  p  (k)e     as  the  operational 

A    J  —  K  ~~ 
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form  of  the  image  source.   Second,  we  will  use  the  mirror  image  condition, 


as  in  (2.17),  to  satisfy  the  boundary  condition  at  x»  =  0.   Then  in  k  -  space, 
from  (2.13),  we  have 

i  k  T  -  i  k  T 

(elkl2  +  £2k22  +  €3k32)  °R  =  pR(-)e   3   +  pR(kl'  k2'  "k3)e  '       (2*21) 


Let  us  define 


ks  ■  VelV  +  ^2dl^l  (2-22) 


Then  because  of  the  image  condition  the  k  -  plane  singularities  of  interest 

are  at  k  =  +jk  .  As  before  we  evaluate  an  ordinary  Fourier  integral  to 
3      s 

obtain  the  mixed  -  space  representation  of  $  in  Region  2: 

R 

-k  T 

PR(kV  k?)e 

Vkl'  V  x3}  = FT^ Sinh  ksx3  (2'23) 

03s 

where  p  (k  ,  k  )  =  p  (k  ,  k  ,  k  =  jk  ) .   It  is  easily  shown  that  the  0 

i\±idri±djS  R 

given  by  (2.23)  satisfies  the  boundary  conditions  at  x  =0. 

Now  we  apply  the  boundary  conditions.   It  is  clear  that  (2.7)  and 
(2.8)  transform  directly  into  mixed  -  space  representation,  so  that  at 
x  =  -T  we  have 


S3°L(kl'  k2'  "T)  "  e3a30R(kl'  k2'  "T)  =  a(kl'  V^o         (2'24) 


and 


0L(kl9  k2,  x3  =  -T)  =  ^(k^  k2,  x3  =  -T)  (2.25) 
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From  (2.20)  and  (2.23)  we  obtain  the  required  quantities  for  the  mixed  - 
space  boundary  conditions  above,  and,  upon  substitution,  (2.24)  and  (2.25) 
reduce  to 


and 


-k  T       -k  T 
pTe     +  2pDe  s  cosh  k  T  =  2a 

L  K  S 


-k  T         -k  T 
e3ks  pLe     =  2krPR6     Sinh  ksT 


From  these  two  equations  we  obtain  expressions  for  p  and  p  which  give  the 

R      L 

desired  solutions : 

k  (x  +  T) 
c(k   k  )   e  r  5  sinh  k  T 

$L(*l9   k2,  x  )  =     -      k  sinh  k  T  +  e  k  cosh  k  T       (2'26) 

o      r      s     3  s      s 

a(kl5  k2)  sinh  kx 

$R(k1,  k2,  x  )  =       -      k  sinh  R  T        cQsh  k  T    (2.27) 

o      r      s      3  s      s 

Since  the  potential  $  expressed  as  (2.26)  in  Region  1  and  (2.27) 
in  Region  2  satisfies  the  appropriate  boundary  conditions  and  differential 
equations,  by  the  uniqueness  theorem  this  is  the  solution  for  the  static 
fields  in  the  light  -  valve  configuration. 

We  will  not  attempt  to  evaluate  the  x  -  space  representation  of  $ 
except  for  the  simple  charge  pattern  consisting  of  a  constant,  say 
o(x1,x2)  =  A.   Then 


o(k1?  k2)  =  (2n)2  6(k1)  6(k2)A 
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so  that  we  may  obtain  $(x)  by  determining  the  limit  of  0(k  ,  k  ,  x  )  as 
k  ,  k  vanish.   If  this  is  done,  we  find  that 


V-}  =  "X3  A/£o€3  (2*28) 


and 


0L(x)  =  TA/eo£3  (2.29) 


As  might  be  expected,  this  potential  is  identical  to  that  of  an  infinite 

parallel  -  plate  capacitor  filled  with  a  material  of  relative  dielectric 

constant  e  . 

To  conclude  this  section  we  will  comment  on  some  of  the  immediate 

implications  of  the  potential  $  (k  ,  k  ,  x  )  expressed  by  (2.27).  Recall 

R   1    2    3 

that  the  relative  frequency  response  of  a  system  is  given  by  the  Fourier 

transform  of  the  system  response  to  a  unit  impulse  function.   The  unit 

impulse  function  for  a  two  -  dimension  spatial  pattern  is  &(x  )   &(x  ),  which 

is  also  a(xn ,  x  )  for  a  unit  point  charge.   For  a  unit  point  charge, 

a(k  ,  k  )  =  1,  so  $  (k  ,  k  ,  x  )  with  a  =  1  is  the  relative  spatial  frequency 
J.   c.  K  1   2   3 

response  of  the  potential  function  in  any  plane  x  =  constant  (in  Region  2) . 
For  example,  at  x0  =  -T  we  have 


3 


sinh  k  T 
s 


o^v  1'   2'   3       "  k  sinh  k  T  +  eJt  cosh  k  T 

r       s     3  s       s 


This  is  a  two  -  dimension  spatial  frequency  response.   The  frequency 
response  to  a  pattern  of  lines  is  given  by  <J>(k  ,  k  )  evaluated  along  a  line 
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in  the  k  ,  k  plane.   If,  in  general,  e  ^  e  then  the  value  of 

k  =  Ve  k  d        :   k  ^/TiT  (and,  in  turn,  the  relative  line  response)  will 

S       J.  X       d.    d.  j 

depend  on  the  direction  of  the  k  ,  k  line  (which  is  normal  to  the  direction 

of  the  x  ,  x  lines).  However,  if  we  have  a  uniaxial  crystal  (such  as  KDP) , 

then  £     =  £     and  the  frequency  response  to  a  line  pattern  is  given  by 
(setting  k  =  0,  say) 


*(kr  0,  -T)  =  -± 


sinh  k  T1 


e  kn  sinh  k  T'  +  VeTT,,  cosh  k  T' 
o  1       1      13      1 


where  T'  =  T/e, /e  .   If  this  $(k  )  is  normalized  (i.e.,  $(k  )/$(k  =  0))  and 
the  hyperbolic  functions  changed  into  exponentials,  the  resulting  expression 
will  be  very  similar  to  the  function  V(N)/V(0)  discussed  in  Section  I,  if  we 
identify  k  I1  in  $(k  )/$(0)  with  the  exponent  jtNT  in  V(N)/V(0).  The  critical 
difference  is  that  the  equivalent  crystal  thickness  I"  =  TveT/eT  is  a 
function  of  the  ratio  of  the  relative  dielectric  constants  e-.  and  e^.  For 


,C3] 

-•1  ~~  ^3 


example,  using  published  values    of  €-,  and  eo  for  some  typical  light  - 
valve  crystals  we  have 


KH  P0,   (KDP)  T'  ~  l.UT 

2  4 

KD  P0,   (KD*P)  T'  ~  1.1T 

NH,  H  P0,   (ADP)         T*  ~  1.9T 

Thus,  for  a  given  crystal  thickness,  the  normalized  frequency  response  of 
the  potential  function  should  be  best  for  KD*P,  for  which  the  effective 
thickness  is  nearly  equal  to  the  real  thickness.   However,  this  does  not 
imply  that  of  the  three  crystals  KD*P  is  best  for  use  in  a  light  valve,  since 
we  have  not  considered  the  comparative  effect  of  the  electro-optic  co- 
efficients of  each  crystal. 
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III.   FIRST  -  ORDER  OPTICAL  SOLUTION   • 

Materials  such  as  crystals  which  are  electrically  anisotropic 
for  static  and  low  -  frequency  fields  are  also  anisotropic  for  optical  - 
frequency  fields.   Thus  the  optical  electric  material  equation  is  also  a 
tensor  relation,  although  the  optical  dielectric  constants  may  differ 

considerably  from  the  low  -  frequency  dielectric  constants,  since  dielectric 

fpl 

properties  are  frequency  -  dependent.   Further,  it  can  he  shown     that  the 

electro-optic  effect  can  be  treated  by  introducing  additional  terms  into  the 
optical  dielectric  tensor.   These  electro-optic  terms  are  small  compared  with 
the  "natural"  dielectric  coefficients,  so  that  they  may  be  treated  as  a 
perturbation  effect.  Using  this  approach  we  will  obtain  a  first  -  order 
perturbation  solution  for  the  propagation  of  light  through  the  electro  - 
optic  crystal  in  the  light  valve. 

The  linear  electric  material  equation  is  given  by  the  tensor 
relation 


D  =  eo  e  •  E  (3-1) 


where  e   is  the  relative  dielectric  tensor.   In  the  absence  of  the  electro  - 
optic  effect  the  elements  of  e  are  constants  (the  "natural"  or  unperturbed 
dielectric  coefficients).   To  express  the  electro  -  optic  effect  additional 
terms  (functions  of  the  static  field  present)  are  added  to  the  natural  di- 
electric coefficients.   The  functional  expression  of  these  electro  -  optical 

r  21 

terms  depends  on  the  symmetry  class  of  the  crystal  being  considered.     We 

will  consider  only  the  symmetry  class  ^+2m  (which  includes  KDP  and  its 
isomorphs),  for  which  e  has  the  following  form  in  the  principal  -axes 
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coordinate  syst 


(3.2) 


The  constants  n  and  n  are  respectively  the  ordinary  and  extraordinary 
indices  of  refraction  (n,  =  n  so  that  ^+2m  crystals  are  uniaxial) .   The 
quantities  a.  are  the  linear  electro  -  optic  terms  and  are  related  to  the 
static  field  E  and  the  linear  electro  -  optic  constants  r.  ,  ,  r^  by-2-! 


22 
ai  =  -ni  n3  r4lEsl 

a2  =  -nx2  n32  rUl  Es2  (3-3) 

k 

3      1   63  s3 


In  the  light  valve  the  static  field  is  in  general  a  function  of  position  so 
that  the  a.  in  (3-3)  are  not  simple  constants.   In  particular,  the  static 
field  is  determined  by  the  potential  0  given  by  (2.27): 

$(1^,  k2,  x  )  =  p(kl5  kg)  sinh  kgx  (3-*0 


for  -T  <  x_  <  0  and  where  we  have  defined  p  as 

-  a(k   ,   k   )/e 

aflr  k     )     II     L  d  O (^.5) 

P^j    ^2)    -     k     sinh  k  T  +       k     cosh  k  T  VJ  •; 

r  s  3  s  s 
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By  definition  E  =  -  v§>  so  that  in  mixed  -  space  representation  (where  the 


— s 


operator  o\  is  unchanged  but  o\  -*  jk  ,  d0  -»  jk  )  we  have 


2  2 

a1(k1,  k2,  x  )  =  jk^  n3  r^  p  sinh  k^ 

2   2 
a2(k1?  kg,  x  )  =  jk^  n3  r^  p  sinh  kgx  (3-6) 

a3(k1,  k2,  x3)  =  k^  r63  p  cosh  k^ 


As  mentioned  Section  II,  we  will  not  need  the  x  -  space  representation  of 
the  functions  a.  above,  since  the  analysis  to  follow  uses  the  a.  in  mixed  - 
and  k  -  space  only. 

It  is  convenient  to  separate  the  dielectric  tensor  e  into  two 
terms,  one  containing  the  constant  coefficients  n.  and  the  other  containing 


the  functions  a. ,  so  that  (3-2)  becomes 


|  =  |n  +  l  (3-7) 


where 


In  =[0    n/  0   |  (3.8; 


and 
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0    a3  a 

1  -Ia3   °   a] 
a2   &1  0 


L 


(3.9) 


To  obtain  the  appropriate  electromagnetic  field  equations  we  must 
start  with  Maxwell's  equations.  We  assume  all  field  and  source  quantities 
are  time  -  harmonic,  so  that,  for  example,  the  electric  field  vector  E  is 
given  by 

E(x,  t)  =  E(x)e_Ja3t  (3.10) 

Then  Maxwell's  equations,  relating  the  usual  field  vectors  E(x) ,'  D(x) , 
H(x),  and  B(x)  to  the  source  current  density  J(x)  and  the  source  charge 
density  p(x),  are  given  by 

V  x  H  =  J  -joD  (3-11) 

V  X  E  =  jouB  (3-12) 

V  •  B  =  0  (3-13) 

V  •  D  =  p  (3-1*0 

By  combining  the  divergence  of  (3. 11)  with  (3.1*0  we  obtain  another  useful 
equation,   the  equation  of  source  continuity: 

V  •  J  =  joop  (3.15) 
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In  addition  to  Maxwell's  equations,  we  must  use  the  material 
equations  to  relate  D  to  E  and  B  to  H.   The  electric  material  equation  is 
given  by  (3-1)?  and  the  magnetic  material  equation  for  most  crystals  (in 
particular,  KDP  and  its  isomorphs)  is  simply  the  free  -  space  permeability 
equation 


B  =  h-oS  (3.16) 


The  above  equations  together  with  appropriate  boundary  conditions 
are  sufficient  to  determine  any  of  the  field  vectors.  Further,  if  any  one 
field  is  known,  say  E,  the  other  three  are  given  directly  by  the  first  two 
of  Maxwell's  equations  and  the  two  material  equations.   This  known  field  is 
the  solution  of  a  wave  equation  that  relates  the  single  vector  to  the 
sources.   In  our  case  the  most  convenient  wave  equation  involves  E.   One 
such  wave  equation  for  E  is  obtained  by  an  appropriate  combination  of  (3. 11), 
the  curl  of  (3. 12),  and  the  material  equations  (3*1)  and  (3-l6): 


V  X  V  X  E  -k  '"  1  •  E  =  J041  J  (3.17) 


2    2 
where  k   =  co  \i  e    .      However,  this  wave  equation  is  difficult  to  analyze, 

because  of  the  term  V  x  V  x  E.   By  use  of  the  vector  identity 


V  X  V  X  E  he  V(V  •  E)  -  S^E 


we  obtain 
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V  E  +  k   €  •  E  =  -joun  J  +  V(V  •  E) 

an  appropriate  combination  of  (3-1) ,    (3.1*0 3  (3-15)  and  the  expression 
(3.7)  f°r  §  we  have 

\       *  ■   £  ■  l^f     -  V  •  1  •  £  -  (n32  -  n±S)   93E3 

When  this  expression  for  V  •  E  is  substituted  into  the  wave  equation  we 
obtain 

2 

^  +  ko2ln  •  e  +  (^  - 1)  va3E3 

=  "[4  +        2     2     W]    '    [(jC4io-  +  ko\    '   -]  (3'l8) 

n.    k 
1     o 

where  I  is  the  identity  tensor.   For  our  purposes  the  wave  equation  (3 .18) 

is  more  convenient  than  (3.17).   Also,  in  free  space  n  =  n  =  n  =  1  and 

£  vanishes,  so  that  the  free  -  space  wave  equation  corresponding  to  (3.18)  is 

^  E  +  kQ2  E  =  -jo^o  [I  +  -^  W]  •  J  =  M,  say,  (3-19) 

k 
o 

where  the  "modified"  source  current  M  is  defined  to  be 

M  =  -jcm  [I  +  -^=  W]  -J  (3-20) 

k 
o 

Equation  (3.19)  represents  an  uncoupled  partial  differential 
nation.   That  is,  (3 .19)  is  three  separate  wave  equations,  one  for  each 
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component  of  E.   However,  the  crystal  wave  equation  (3.18)  is  coupled 

because  of  the  terms  involving  ^  •  E  and  VcLE  .   Further,  since  the 

coefficients  of  y   are  functions  of  position.  (3.18)  is  a  linear  differential 

equation  with  non  -  constant  coefficients,  so  that  the  solutions  must  be 

nonlinear  functionals  of  the  elements  of  £.   In  the  light  valve  this  means 

that  the  visible  picture  is  a  nonlinear  function  of  the  applied  charge 

pattern.  As  will  be  seen,  however,  the  first  -  order  perturbation  solution 

is  linear  in  form. 

We  have  stated  that  the  elements  of  the  effective  electro  -  optic 

tensor  y   are  small  relative  to  those  of  the  dielectric  tensor  €  •   This  is 
^  =n 

true  if  the  magnitude  of  the  static  field  is  reasonable.   For  example,  a 
typical  light  -  valve  crystal  has  e  ~  25,  n  ~  n  ~  1.5,  and  r^-  ~  r,   ~ 
10    m/y.   Then  for  a  typical  crystal  surface  potential  ~5kV  and  crystal 

thickness  ~0.2mm,  the  magnitude  of  the  static  electric  field  is  ~10  v/m,  so 

-5 
that  an  element  of  y   is  ~10   relative  to  an  element  of  e  .   Thus  x  is  indeed 

^  =n       = 

small  compared  with  e   ,    and  a  first  order  perturbation  solution  should  give 

a  reasonable  approximation  to  the  exact  solution. 

T151 
Recall  that  the  perturbation  approach     may  be  applied  to  a 

two  -  operator  equation  of  the  form 


H  •  E  -   N  +  \G   -   E 

which  relates  the  solution  vector  E  to  the  source  vector  N.   The  operator 
H  represents  the  "unperturbed"  system  (for  which,  presumably,  all  solutions 
are  known  or  can  be  obtained),  and  the  operator  G  represents  the  system 
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perturbation.   The  parameter  X   is  a  small  constant  representing  the  "strength" 
of  the  perturbation  (the  perturbation  is  "turned  on"  as  \   is  increased  from 
0  to  1).   The  first  order  perturbation  solution  for  X  =  1  is  then  given  by 


E  ~  So  +  h 


where  EL  is  the  solution  for  the  unperturbed  system 


and  E,  satisfies 


H.^-H 


H  •  E_x  =  G  •  EL 


We  will  refer  to  EL  as  the  zero  -  order  term  and  E,  as  the  first  -  order 
term  of  the  complete  first  -  order  solution.   Note  that  the  source  N  and 
the  uniqueness  theorem  preclude  the  possibility  that  the  unperturbed  solution 
EL  is  degenerate. 

In  the  absence  of  the  electro  -  optic  effect  (when  the  static 
field  vanishes  so  that  ^  =  0)  the  crystal  wave  equation  (3«l8)  becomes 

2 

^E  +  k  2  en  '  E  +  (-2-  -  1)  ?c)E 
-    o  =n   -      2       3  3 

=  -W0  [I  +  ~Y~2     W]  '  ^  (3-21) 

n  k^ 
1  0 

=  N,  say, 
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where  the  modified  crystal  source  current  N  is  defined  to  be 


N  =  -ja*i  [I  +  -~2   W]  *  ^  (3'22) 

n.  k 
o 

Evidently,  the  unperturbed  system  operator  H  is  the  partial  differential 
operator  on  the  left  of  (3.21),  so  that  the  zero  -  order  term  E  satisfies 


(3.21).  Further,  the  perturbation  operator  G  is  clearly 

-[k  2I  +  -\     W]  •  x 

nl 

so  that  the  first  -  order  term  E  satisfies  (3.21)  if  the  "zero  -  order 
source  current"  jcou  J  is  replaced  by  the  "first  -  order  source  current" 

The  uniaxial  wave  equation  (3.21) only  partially  describes  the 
unperturbed  system,  since  it  governs  only  the  fields  inside  the  crystal.  All 
fields  and  sources  outside  the  crystal  must  satisfy  the  free  -  space  wave 
equation  (3.19)  ■>    and  the  dielectric  discontinuity  at  the  crystal  surfaces 
introduces  boundary  conditions.   The  optical  configuration  of  the  crystal 
(see  Figure  3)  is  similar  to  that  for  the  static  field:   the  geometry  is 
identical  but  the  right  -  hand  interface  is  not  a  conducting  plane  (the  con- 
ducting film  --such  as  vacuum  -  deposited  gold —  is  essentially  transparent 
to  light  ) .   We  will  use  the  same  coordinate  system  as  before  and  again  assume 
the  crystal  extends  to  infinity  in  the  x  x  plane.  Regions  1  and  2  are 
defined  as  before,  and  Region  3  denotes  the  right  half  -  space  x  >  0.  How- 
ever, we  use  subscripts  "L"  and  "R"  to  denote  variables  for  Regions  1  and  3, 
respectively,  and  the  subscript  "C"  for  Region  2  variables. 
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%' 

Figure  3.   Crystal  Configuration  for  Optical  Fields. 
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Accordingly,  the  boundary  conditions  are  as  follows:   continuity 
of  the  normal  component  of  D  and  the  transverse  component  of  E  requires 


n3  EC3  =  Ef3 


EC1  =  Efl'  EC2  =  Ef2 


3.23) 


at  x  =  0  or  x  =  -T  (with  the  subscript  "f"  taken  as  "R"  or  "L").   Con- 
tinuity of  H  at  the  interfaces  requires 


V  X  E^  =  V  X  Ef  (3-24) 


at  each  interface. 


The  Zero  -  Order  Term.   The  zero  -  order  term  E  can  easily  be 
obtained  by  conventional  methods,  since  the  "source"  is  an  x  -  directed 
plane  wave  in  Region  1,  so  we  have  normal  incidence  at  the  crystal  interface 
x  =  -T.   Recall  that  for  an  isotropic  dielectric  sheet  we  assume  that  the 


solution  vectors  for  a  plane  wave  at  normal  incidence  have  no  x-3  -  components 
Assuming  then,  that  E  _  has  no  j 
wave  equation  (3.21)  reduces  to 


Assuming  then,  that  E__  has  no  x  -  component,  we  find  that  the  uniaxial 

UL         3 


VEM  +  n,  \ 


"1  "o  -0C   ° 


where  N  =  0  since  the  source  is  confined  to  Region  1.   This  is  identical  to 
the  isotropic  wave  equation  so  that  the  analysis  and  solution  for  the 
isotropic  plane  sheet  are  valid  here.   This  solution1-    is  indicated  in 


3^ 


Figure  k   and  is  as  follows .   In  Region  1  we  have  the  incident  plane  wave 

E.    =  E  e   '3,  where  E  is  a  constant  vector  and  E   •  x„  =  0.   The 
-mc   -a  -a  -a   -3 

reflected  wave  E   „,  is  of  the  form  e     3  and  is  related  to  the  incident 
— refl 

wave  by  a  (reflection)  coefficient  r,  ,  where 

jL 

,   2   ,v  n    -j2nnk  Tn 
iPk  T    (nl   -  1)  [1  "  e  d  1  o  ] 

r  =  e  3dKoL   ± (3.25) 

L  /      ,  ,2  -j2n,k  T   ,     nx2  K3"&J 

(n1  +  1  )      e    1  o  -  {n±   -  1) 

so  that 


Sol  "  5a  C^*3  +  rL  e"Jk°X3l  <3.26) 


The  transmitted  wave  E.     is  of  the  form  e    3  in  Region  3  and  is 

— tran 

determined  by  the  (transmission)  coefficient  rD, 

Fr  =  („x  +  I)2  e"^1  -  (ni  -  I)2  ejnlkoT  (3,27> 

so  that 


£or  "  Sa  rR  edk°X3  (3.28) 


In  the  crystal,  Region  2,  we  have  a  positive  -  directed  and  a  negative  - 

directed  plane  wave,  with  wave  number  n.k  .   The  solution  for  E--  is  then 

l  o  — uu 

given  by 

^0C  =  EQ^c[eJnlkox3  +  T!e-dnlkox3]  (3.29) 
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% 

Figure  4.   Zero  -  Order  Optical  Solution. 
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where 


(3.30) 


and 

2(n  +  1)  e"JkoT 
"C  =  (nx  +  I)2  e-^lkoT  -  (n]_  -  if   e^lV  (3>31) 

It  is  not  difficult  to  show  that  the  solution  for  the  zero  -  order  term 
given  above  satisfies  all  the  conditions  required. 

The  First  -  Order  Term.   The  solution  for  the  zero  -  order  term 
was  easily  obtained  with  an  assumption  suggested  by  the  equivalent  isotropic 
configuration.  No   such  analogy  is  evident  when  the  source  is  inside  the 
crystal,  as  in  the  case  of  the  first  -  order  term.   However,  a  solution  can 
be  found  by  a  more  general  approach.   For  example,  the  method  of  images  is 
appropriate  here  because  of  the  form  of  the  boundary  conditions. 

We  will  first  obtain  the  general  form  of  the  solution  in  the  free  ■ 
space  Regions  1  and  3-   In  k  -  space  the  wave  equation  (3. 19)  is 


(k2  -  kQ2)  E(k)  =  -  M(k)  (3.32) 


is  an  image  source,  the  mixed  -  space  representation  of  E  is  determined 
by  the  singularity  at  k  =  +k  ,  where 

kf  =  (ko2  -  kr2)l/2'   kr2  =  kl2  +  k22  (3'33) 
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Since  both  k  singularities  are  on  the  path  of  integration,  we  invoke  the 

radiation  condition  (which  requires  that  a  source  emit  waves)  to  select 

the  desired  singularity.  Accordingly,  k  =  -k  represents  waves  propagating 

in  the  -x0  direction,  so  that  the  general  solution  for  E  in  Region  1  is 
~3  h 

given  by 


-jkfXo 
EL(k1,  k2,  x3)  -  ^  -—-  (3.3k) 


f 


and  in  Region  3  we  have 


Vkl'  k2'X3}  =^W  (3-35: 


f 


where 


«l  "  4(ki'  k2>  k3  =  "V 
4  =  Vkl>  k2>  k3  "  V 


!3.36) 


We  are  ultimately  interested  in  an  expression  for  E  ,  the  first  -  order  term 
for  the  transmitted  wave. 

The  general  forms  of  E  and  E  and  the  boundary  conditions  can  be 

— R      L 

used  to  obtain  a  set  of  conditions  which  the  crystal  field  EL,  must  satisfy 

— u 

at  the  interfaces.   Since  there  are  no  sources  in  Regions  1  and  3?  the  free  - 
space  fields  must  satisfy  V  •  E  =  0  everywhere  in  the  appropriate  regions. 
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Thus,  evaluation  of  7  •  E  =  0  (in  mixed  space)  at  the  interfaces,  followed  by- 
appropriate  substitutions  from  the  boundary  conditions  (3-23)  yields 


klECl  +  k2EC2  "  n32  kfEC3  =  °  at  X3  =  'T 

(3-37) 
klECl  +  k2EC2  +  n3  Vc3  =0atV° 


If  we  evaluate  the  free  -  space  components  of  the  boundary  condition  (3.24) 
and  make  appropriate  substitutions  from  (3-23)  we  obtain  four  more  indepen- 
dent relations  for  the  interface  values  of  E^,  of  which  two  contain  E   : 


jkx(n32  -  1)EC3  =  -(jkf  +  B3)EC1  at  x3  =  -T 

(3-38) 
jkx(n32  -  1)EC3  =  (jkf  -  ^3)EC1  at  x3  =  0 


The  other  two  equations  are  similar  and  relate  E   and  E   .   An  appropriate 
combination  of  these  equations  (3-38)  yields 


(kf  -  jd3)(k2Ecl  -  k±Ec2)    =0  at  x3  =  -T 

(3-39) 
(kf  +  j£>3)(k2Ecl  -  k^)  =0  at  x3  =  0 

Another  combination,  of  the  equations  in  (3«38)  and  those  in  (3.37)?  results 
in 


(n32kQ2  -  k/  -  jn32kfa3)(klEci  +  k2EC2)  =0  at  x3  =  -T 

(3^0) 
(n32kQ2  -  kr2  +  jn32kfa3)(k1Ecl  +  k^)  =0  at  x3  =  0 


To  obtain  the  crystal  field  EL,  we  need  the  k  -  space  uniaxial 
wave  equation  which,  from  (3.21),  is 

2 

k2^  -  k  2e   •  E_  +  (-^-  -  1)  k  k  E 
-C    o  =n   -C    K      2  ~  3  C3 

n,  J 

1  (3.^1) 

^o       2   2 

=  "I  =   2  2  Cni   o  -  _  -(^  '  ^)] 
nn  k 
1  o 

It  is  possible  to  manipulate  (3. 1+1)  such  that  three  equations  are  obtained, 
one  for  each  component  of  EL.  With  these  equations  one  may  easily  obtain 
two  useful  expressions: 


Al  S  k2ECl  "  R1EC2  =    t       A  0-te) 

k3  "ko 


U         k  •  K 

A2  =  klECl  +  k2EC2  "  7  2  \  2  --2—2        (3'U3) 

k3  "  ke     k3  "  kc 


where  for  convenience  we  have  defined 


/  2.  2   .  2x1/2 

k  =  (n  k   -  k  )  ' 
c     1  o     r   ' 

k.  -  <Vty<»3V  -  kr2>1/2 


(3. mo 


The  definitions   of  A     and  A     are  appropriate  because  of  the   form  of  the 


1+0 

interface  conditions  (3«39)  and  (3-^0)  and  because  these  expressions  are 

useful  in  the  analysis  to  follow. 

The  image  source  solution  requires  two  image  sources,  one  for  each 

of  the  image  regions  x  <  -T  and  x  >  0.   Until  the  boundary  conditions  are 

applied,  the  two  image  sources  must  be  treated  as  independent  functions  and 

have  the  following  characteristics:   the  left  -  hand  image  source  N  vanishes 

for  x  >  -T  and  produces  positive  -  directed  waves  in  Region  2,  and  the 

right  -  hand  image  source  N  vanishes  for  x  <  0  and  produces  negative  - 

- K  3 

directed  waves  in  Region  2.   These  conditions  on  the  image  sources  are  required 
since  Maxwell's  equations  are  linear,  and,  in  general,  conditions  at  one 
boundary  may  be  independent  of  those  at  the  other. 

For  each  of  A  and  A  we  have  A.  =  A^      +  A  ,  where  A  .  is 
determined  by  the  true  source,  and  A   by  the  image  sources.   If  we  define 


kcPL  "  k2NLl  "  klNL2 


and 


kcPR  "  k2NRl  "  kA2 


then  A  ,  obtained  from  (3.1+2)  has  the  mixed  -  space  representation 


2jAI1  =  PLeJkcX3  +  PRe"Jkcx3  (3.I+5) 


Expressions  for  P  and  P  in  terms  of  the  source  function  A   may  be 


in 


obtained  from  the  interface  conditions  (3 • 39) 3  by  substitution  and  evaluation 

of  ATn   as  given  above.  However,  we  only  need  A,...  at  x  =0,  for  which  we 
II  II     3 

have,  from  (3.U0)  and  (3-1+5) 

p     +p     _     2.   ^1^    +  Cl(°^kc   -  kf)eJk°T   +  K  +  V6'^        (3.46) 
L         R  (k     +  kj2  e"JkcT  -    (k     -  kj2  ejkcT 

where 


C-^-T)   =  (kf  -   jd   )   A^     at  x3 


=  -T 


and 


^(0)   s   (kf  +  j^3)   AN1     at  x     =  0 


The   function  A,    =  A       +  A       is,   then,    at  x     =  0: 


Al  ■  \L  +  5J(PL  +  PR)  (3.W 


where  (P  +  P_.)  is  obtained  as  above. 

L     K 

To  determine  the  appropriate  expression  for  P^      from  (3-1+3) ,  we 
note  that  by  definition  in  (3. 1+1),  the  quantity  k  •  N  reduces  to 


-JwM-0      o  2    2s 


k  •  N  =  — g  (n1  kQ"  -  k")(k  •  J) 


n_,  k 
1  o 


1+2 


so  that  an  alternate  expression  for  (3.^+3)  is 

k  N        jujx 

A2  =727^--2-2-  '  "  (3.W) 

k3  "  ke    nl  ko 

For  the  image  function  A-r9,   we  require  that  the  image  sources  vanish  in 

Region  2,  so  the  term  in  k  •  J  does  not  affect  the  mixed  -  space  representation 

of  A.,.-  in  Region  2.   Accordingly,  we  have 

-2JAIp=NL3e    3.^3   e3 
so  that  evaluation  of  the  interface  conditions  (3*^0)  yields 


NL3  ■  NR3  (3'U9) 

2n12kfC2(-T)    +  C    (0)[(k     +  ni2kf)e       e     -    (k     -  n  \   )e     e    ] 
=  2j  


(ke  +  n1\)2e-^T  -    (ke  -  n;L2kf)2eJke- 


where 


C2(-T)  s  r  K2  -  J-^A' V  at 


x„   =  -T 


and 


C2(0)  =  i  (ke2  +  jn12kfB3)AN2  at  x3  =  0 
e 


Then  we  have  Ap  =  A^p  +  A   evaluated  at  x  =  0: 


1*3 

A2  "  V  "  Ti   <NL3  -  V  (3'50) 


where  (NT 0  -  N  )  is  given  above. 


To  evaluate  the  iL  .  we  need  the  first  -  order  current  source, 
determined  by  the  dot  -  product  of  the  electro  -  optic  tensor  (3«9)  and  the 
zero  -  order  crystal  field  (3-29).   With  the  following  definitions 


si  s  rcpni%%i 
g3  s  r<P\r63 


(3.5D 


and 


1    a2 


F2  s  Eal  (3-52: 


F  =  k  F  +  k  F 
3    11    2  2 


and  using  (3 .6)  we  obtain  the  x  -  component 


jn  k  x      -jn  k  x 
(&  •  EQC)1  =  Fxg3ks  cosh  ksx3(e  x  °  -J  +  ne   X  °  ^)        (3-53) 


The  x  -  component  is  given  by  a  similar  expression,  with  F  replacing  F  . 
For  the  x  -  component  we  have 


hk 


jn  k  x      -jn  k  x 
S  •  E^y  J*3&±   sinh  ksx3(e  X   °  3  +  ne   X  °  3)         (3-5*0 


;e  the  electro  -  optic  effect  is  confined  to  Region  2  and  vanishes 
elsewhere,  the  k  -  space  representation  of  ^  •  IjL_  is  obtained  by  a  x 
Fourier  integral,  integrating  from  x  =  -T  to  x  =  0. 

Such  a  "truncated"  Fourier  integral  is  equivalent  to  a  convolution 
integral  of  the  form 

-ikx    P°°   dk'   eJT(k  "  k'}   1 


-00 


where  f (k ' i  is  the  ordinary  Fourier  integral  (2.11)  in  one  dimension  only. 
We  will  denote  the  above  convolution  integral  by 

R(k  -  k')  *  f(k') 


where 


R(k)  =  (eJkT  -  l)/jk 


With  the  definition 


loo      -jk  x 


poo       -JK  x 

L(k)  s  /   dx3e   *   3  £(kx,  k2,  x3)  '  EQC(x3) 

—  00 


the  first  -  order  source  current  becomes 
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ja*i  J(k)  =  k  2R(k  -  k»  )  *  L(k')  (3-55) 

o o    3     3 


where  k'  =   (k  ,  k  ,  k'  ).  Then  from  (3-53)  we  have 

Lx(k)  =  «Flg3ks[8(k3  -  k*)  +  6(k3  -  ka)  (3.56) 

+  n6(k3  +  ka)  +  n6(k3  +  k*)] 
where 


k  =   n_,k  +  j'k 
a    1  o     s 


A  similar  expression  is  obtained  for  L  ,  and  from  (3-54)  we  have 


L3(k)  =  gnF3gl[5(k3  -  ka)  -  5(k3  -  ka)  +  n&(k3  +  kQ)  -  n6(k3  +  k*)  ]    (3-57) 

We  have  defined  the  first  -  order  source  current  as  the  convolution  (3. 55) 
because  the  delta  functions  in  L  reduce  the  convolution  to  a  simple  sum  over 
the  appropriate  values  of  k'.   If  a  mixed  -  space  representation  is  to  be 
obtained,  the  k'  convolution  may  be  deferred  until  after  the  k  transform 
of  R(k  -  k')  is  obtained. 

We  now  evaluate  the  mixed  -  space  expression  for  A  .  From  (3.4-1), 
(3-42),  and  (3-55)  we  have 

R(k  -  kf) 

ANl(-)  =   o  1 2  *  ^iV-')  "  k2Ll(-')]  (3"58) 

k3  "  kc 


U6 


If  we  define 


Gc(k,  k')  =  -R(k3  -  k^)/(k32  -  kc2) 


and  evaluate  the  mixed  -  space  representation  at  the  interfaces  we  obtain 

JT(kc  -  k') 

G  (k'  ,  x  -  0)  =  ^Z-2 

3        2k  (k  -  k') 


cv  c    3' 

-jk'T    jkj 

i te',  x  =  -T)  =  ^ =-S 

^  2k  (k  +  kl) 

c  c    3' 


and  the  derivatives  reduce  to 


a3GQ(k',  x3  =  0)  =  JkcGc(k',  x3  =  0) 


a3Gc(k',  x3  =  -T)  =  -JkcGc(k',  x3  =   -T) 


The  relation  between  the  interface  value  of  cVG  and  that  of  G  is  independent 

3  c  c 

of  k',  so  that  after  convolution,  the  same  relation  holds  between  A   and 
^X,  at  the  interfaces.   This  result  reduces  the  expressions  (3.^6)  and 
(3.V7)  to 


Ax(k  ,  k  ,  x  =  0)  = 

-  ik  T 

2k  (ke  -  kf>V-T> +  'kc  +  VV°>e    c       (3-59) 

"  c  (k     +  kf)2e",'kcT  -   (k0  -  kf)2edkcT 


^7 


where 


V-T> E  Vki>>vx3  ■  "T) 


with  a  similar  expression  for  A^, (0),  at  x  =  0. 

The  exact  expression  for  A  at  x  =  0  is  then  obtained  by- 
substituting  into  ( 3- 59)  the  interface  values  of  G  for  those  of  A^ 

and  evaluating  the  k0 '  convolution  as  indicated  by  the  definition  of 

3 

G  and  (3-58) •   The  result  is  the  following  expression: 


A1(k1,k2,x^  =  0)  = 


g3ksko2(klVk2Fl> 


2      2  2  2   "JkcT  2      ^kcT 

V*r^\   +  (r-i)kr  ][(Vkf)  e  -  (VV  e 


-  ik  T 
x    {[kc(7-l)(l+n)   +  n1ko(7+l)(n-i)](kc+kf)e  ' 

ik  T 
-    [kc(7-l)(l+1)   -  niko(7+l)(ri-l)](kc-kf)e     c 

-jn  kQT  Jn  kQT 

-   2k  coshk  T[k   (7-l)(e  +  ne  ) 


-jnnk  T  jnnk  Txn 

i/      ,\/  1  o  °   1  o    ) ] 

"  n1k0^+1)(e  "  ne 


jnxkoT  jn^T 


-   2o(7/ks)sinhksT[2nikokf(e  -  rje  ) 

(0n   2V   2       .2,      _.;      "JnlV  JnlkoT)]} 

+   (2n1  kQ     +  kr    (7-l))(e  +  ne  J2S 


1+8 

where  y   =  £-,/€-,.   As  can  be  seen,  this  exact  expression  is  quite 
complicated.   However,  a  reasonable  approximation  results  in  considerable 
simplification.  We  are  discussing  an  optical  solution,  so  that  the 
free-space  wave-number  k  =  2n/\     is  very  large.   Thus  for  light  of 


o      ^  '         o 


~  0.6um,  we  have  kn  ~  10  m   .   This  value  is  to  be  compared  with 


typical  values  of  k  and  k  .   From  (3. 51)  and  (3*5)  we  find  that  the 
above  expression  for  A, (0)  contains  the  factor  a(k,,k^)  --  the  Fourier 
compos i .     f  the  crystal  charge  pattern  --  so  that  the  meaningful 
range  of  k,  and  k^.  is  determined  by  the  maximum  frequency  content  of 
s(k, ,k  ).  For  standard  television  we  have  a  maximum  resolution  of 

500  lines  (half -cycles)  in  any  direction,  which,  for  a  typical 

k   -1 
crystal  size  of  25mm  square,  gives  ~  10  m   for  a  maximum  value  of 

2     2   2   2 
k  or  k  .   Then  for,  say,  k   =  k  -k  -k  ,  the  approximation  k  ~  k 

J.       c.  1       O     -L    <—  IO 

is  very  good  for  all  values  of  k, ,k  of  interest.   Also,  from  (3.UU), 
we  have  k  ~  n,k  .   These  approximations  reduce  the  expression  for 
A, (0)  to  only  one  term  of  magnitude  k  relative  to  the  remaining  terms 
Retaining  only  this  k  term,  we  obtain 


A1(k]_,k2,x3  =  0)  =  k2Ecl-klEc2  Z 

k_Fn    -  k-P0     .  2   J2nik0T 

(jg_k  sinhk  T)    2  X      .    1  2     -1  +   H  e.0     .    _  (3-60) 

Vde>3  o  s    '        n   4l  2  j2n1kQT 

1  -  q  e 


where,  as  before,  rj  =  (n  -l)/(n ,-»l). 
To  obtain  an  expression  fc 
using  the  first  -  order  source  current  ( 3 • 55) : 


To  obtain  an  expression  for  A  at  x  =  0,  we  begin  with  (3.48), 


1*9 


k  N 
\2  =   2"  P  "—  R*[k-  L  (k-)] 

k3  "  ke     nl 

By  substitution  for  N  ,  from  (3. 1+1),  this  becomes 

n  A  2(k  L  +  k  L  )  -  n A  2(k  •  L) 
A   =  R*  3  ° Li ?_2 1  r   (3.61) 

n3  (k3  -  ke  ) 

When  this  expression  is  evaluated,  the  term  obtained  from  the 

(k,L,  +  k  L  )  component  of  (3.61)  will  be  of  order  k  relative  to  that 

obtained  from  the  k  •  L  component.   This  is  reasonable  since  the  k  •  L 

term,  when  evaluated,  obtains  a  factor  of  order  k   (from  k  L„)  while 
»  '  o       3  3 

2 
the  term  for  (k, L  +  k  L  )  already  possesses  the  factor  k  ,  as  in  (3.6I), 

Since  the  approximation  for  A,,  (3.6O),  contains  only  the  terms  of 

highest  order  in  k  ,  we  may  then  approximate  A^  by 

^2  ~  K    R*  \\        V  (3.62) 

k-   -  k 
3     e 

Note  that  this  approximation  implies  that  only  r^  modulation  is  retained 

in  this  first-order  evaluation,  since  the  factor  r,   occurs  only  in  the 

term  k  •  L.   The  form  of  (3.62)  above  is  similar  to  that  for  A^  in 

(3.58),  and  the  analysis  to  follow  parallels  that  for  A,.   If  we  define 

a  G  (k,  k  ? )  as  in  the  definition  for  G  above  (but  with  k  substituted 
e  ~~   3  c  e 

for  k  ),  the  interface  expressions  for  G  and  cLG  have  the  same  form 
c"  e      3  e 

as  those  for  G  and  S  G  ,  so  that  equations  (3-^9)  and  (3'50)  defining 
A  reduce  to 
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A2(k1,k2,x3  =  0) 


2nx  kf 


2,  <ke  +  nl\>  °  6  V°)  -  <ke  -  "A*  V'1*   (3.63) 


p   2  ~^keT  2   2  ^keT 

(ke  +  nx  Kf )  e      -  (kg  -  n±   kf)  e 


where  A^  (-T)  =  A^  (k  ,k  x  =  -T)  with  a  similar  definition  for  A^(0) 
at  x_  =  0.     before,  we  substitute  G  for  A   in  (3*63)  and  evaluate 
the  convolution  integral.   The  result  is 


A2(k1,k2,x3  -  0)  a 

g3ni2kfks(ko2  '    vkr2)(k1F1  +  k2F2) 

2  2  2  "^keT  2  2  ^keT 

^7nl  ko     +    ^"v^kr    ^nl  kf  +  ke)    e  "    W  kf  "ke^    e  ^ 

-  ik  T 
X    C[(l+i|)(7-v)    -    (l-n)(7+v)(n1ko/ke)](ke+n12kf)e   ' 

ik  T 
+    [(1+t,)(7-v)   +   (l-n)(7  +  v)(n1ko/ke)](ke-n12kf)e     e 

-n  k  T  jn  k  T 

-  2n1coshksT[n1kf(7-v)(e     X  °     +   ne  )] 

-Jn^k  T  jn^T 

-kQ(7+v)(e  -   tie  )] 

-jn  k  T  jn  k  T 

-  2j(7/ks)sinhksT[2n1\okf(e       1  °     -  r\e     X  °  ) 

o     o  o  "jn  k  T  Jnnk  T 

+    (Sn^2  +  kr2(7-v))(e        l0     +  ne     1  °    )]} 
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2   2 
where  v  =  n  /n  .  As  before,  we  use  the  appropriate  approximations, 

k  ~  n,  k  and  k„  ~  k  ,  to  reduce  the  above  expression  to  only  one  term 
e    1  o      f    o' 

of  relative  order  k  .  Discarding  all  but  this  highest  order  term,  the 

o 

result  is 


A2(krk2,x3  =  0)   =  klEcl  +  k2Ec2  S 


k  F     +  k  F 
(-jg3ko  sinh  ksT)     1  I        Is  (3.64) 


Finally,  by  comparison  of  (3 .37)  and  (3-64)  it  is  evident  that 

_  at  x_ 
c3    3 


to  order  k_  we  may  neglect  the  terms  involving  E_n  at  x0  =  0,  so  that 


our  solution  involves  only  E  ,  and  E  ^. 

cl      c2 


The  approximations  for  E  ,  and  E   at  x  =  0  are  easily 
obtained  by  appropriate  combinations  of  (3.60)  and  (3.6^).  For  example, 


we  have 


Ecl(kl'k2'X3   =  0)    =    (k2Al   +  klA2Vkr2 


-jg_k     sinh  k  T 
3  o s 


2  2JnlkoT 


(n1+l)(l-^e       x  u   ) 


2F  k  k  k  2  -  k  2  2jn  k  T 

X    (      2  \  2  +  F  [-+ r-2-     -  n2e       lo]}      (3.65) 

k  k 

r  r 


For  E       the   expression  is   similar  to    (3.65)   differing  only  by  the  inter- 
change of  F,    and  F^  as  well  as  that  of  k     and  k    . 
12  12 
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We  can  obtain  a  more  compact  expression  by  a  change  in 
coordinate  representation,  to  k  ?k  polar  coordinates.   Let 


so  that 


0R  =  tan"1  k2/k1  (3.66) 


k..  =  k  cose,   and  k^  =  k  sine, 
Irk        2    r     k 


Then  the  factor  in  braces  in  (3*65)  reduces  to 


j2n  k  T 

F,  cos26,  +  F0  sin20,  -  Fn  n  e 
1      k    2      k    1  ' 

Further,  recall  (3.52)  relating  F  and  F  to  E^,  which  by  (3-26)  defines 
the  magnitude  and  polarization  of  the  incident  plane  wave.   Let  0 
define  the  angle  of  polarization,  so  that 


F,  =  E  _  =  E  sine 
l   cc2   a    a 


F^  =  E  ,  =  E  cose 
2   al   a         a 


(3.6?) 


(Note  that  elliptical  polarization  can  be  described  by  a  complex -valued 
angle.    J )  This  definition  for  6    ,   together  with  certain  obvious 
triginometric  identities,  allow  a  further  reduction,  so  that  (3 .65) 
becomes 
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-jg  k     sinh  k  T 
E       S  ^_o s 


cl  0  2jn  k  T 

(nx  +  1)(1  -   n  e        ±  °    ) 

2 in  k  T 

x  Ea[sin(2ek  +  ea)  -  n2e  '   1  °    sinej        (3.68) 


Now,  from  (3«^)  and  (3. 51)  we  have 


gn  sinh  k  T  =  n,  iv^r  $(-T) 
&3      s     1  63  c    ' 


where 

<D(-T)  =  0(k1,k2,x3  =  -T) 

This  reduces  (3-68)  to 

-jk  n  r,.0(-T)r  E^ 
F  fk  k  x  -  0^  -  d  o  1  63  c  a 
Ecllkl'k2'X3  "  °j  "  0  2jnk  T 

^  +  1)(1  -  T)  e   X  °  ) 

2jn  k  T 
x  [sin(29k  +  6a)    -   x)   e     °  sin  0^] 

By  a  similar  procedure  we  obtain 

v     rv  „      M  ~    3VlS>3»("T)rcEa 
Ec2(kl5k2,x3  =  0)  =  g  2Jn  k  T) 

(n1  +  1)(1  -    r)  e 

2  2JnikoT 
X    [cos(201    +  e    )   +   n  e  cos  9    ] 


(3.69) 


(3.70) 
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Finally,  from  (3o7)  it  is  evident  that  the  approximation  for 

=  0  is  not  of  order  k  and  so  may  be  neglected. 

o 

We  have  now  obtained  the  approximation  for  EL,  at  x  =  0,  which, 

together  with  the  boundary  conditions  (3«23)  and  the  expression  (3.35) 

for  ED,  determine  the  first  -  order  term  of  the  transmitted  wave  in 
— K 

Region  3.   The  approximation  for  the  first  -  order  solution  in  Region  3 

is  given  by  the  sum  of  ED  as  in  (3 • 35 )  and  the  zero  -  order  transmitted 

— K 

wave  EL       (3»28).  However,  since  the  resulting  expression  would  be 

rather  complicated  (as  in  (3-70) )  we  will  consider  a  special  case. 

If  the  light  -  valve  is  used  to  obtain  visible  pictures 

(intensity  modulation,  then  the  incident  light  is  polarized,  and  the 

transmitted  light  passed  through  an  analyzer  oriented  perpendicular 

to  the  polarization  of  the  incident  light.   In  particular,  if  the 

incident  light  is  polarized  in  the  x,  -  direction,  say,  we  would  be 

interested  in  the  amplitude  (or  intensity)  of  the  E  component  of  the 

2 

transmitted  wave  relative  to  that  of  the  E  component.   If  we  obtain 
the  ratio  of  EQ2,   given  by  (3-70),  to  (EL^)  at  x  =  0,  given  by  (3-28), 
we  should  have  a  good  estimate  of  the  amplitude  modulation  characteristics 

for  this  case.   For  E..  -  polarization  we  have  0   =0  and  F.,  =  E  _  =  0, 

1  *  a  1   a2    ' 

so  from  (3.27)  through  (3-31), 


e^(x  =0)  =  xlEarR  =  Wc  _^l         (3.71) 

n.,  +  1 
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Then,    from   (3-70)    and   (3-71), 


jn A.  k  cJ>(-T) 

(V        /F  ^  I  -  b^    ° 

v   c2'    0Rly|x  =0  0  2jn_k  T 

3  2(1  -    n2e        X  °    ) 

2  2JnlkoT 
X    [cos   28,    +  t]  e  ] 


(3-72) 


where  $(-T)  =  $(k  ,k  ,x  =  -T)  as  given  by  (2.27).  Now  the  far  -  field 

ri  fii 

diffraction  pattern  can  be  shown     to  be  related  to  the  Fourier  decom- 
position of  the  source  pattern,  which  in  our  case  (k  ~  k  )  is  simply 
the  ratio  given  above  (the  ratio  is  used  since  it  is  independent  of  the 
amplitude  and  phase  of  the  incident  light).   If  we  recall  the  discussion 
in  Section  II  concerning  the  static  field  resolution  at  x  =  -T,  the 
function  <|>(-T)  determines  the  pattern  resolution  characteristics  of  the 
static  field.   From  the  ratio  given  above,  $(-T)  is  also  quite  important 
in  determining  the  pattern  resolution  characteristics  of  the  transmitted 
light.   However,  the  presence  of  the  factor 

2  2JnlkoT 

cos  20,  +  r\   e 

indicates  that  the  pattern  resolution  is  a  function  of  direction  in  the 

pattern,  and  that  the  response  should  be  minimum  along  the  direction 

determined  by  0   =  n/h.      This  directional  dependence  is  not  unreasonable, 

since  we  have  observed  that  the  quantity  k  E   +  k  E   is  associated  with 

the  wave  -  number  k  ,  which  in  turn  is  related  to  the  extraordinary 

wave  number  n^k  ,  and  similarily,  kJ5„,  -  k,E_^  is  related  to  the 
3o'  ■"   2  CI    1  C2 

ordinary  wave  number  n  k  .   The  two  components  then  travel  at  different 
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velocities  inside  the  crystal,  and  reach  the  interface  x  =  0  out  of 

2 
phase.   The  result  is  an  amplitude  which  depends  on  the  values  of  k 

2 

and  k   .   Thus  the  pattern  resolution  characteristics  of  the  light  valve 

are,  to  this  first  -  order  approximation,  closely  related  to  those  of 
the  static  -  field  surface  potential,  but  resolution  is  non-uniform  as 
a  function  of  k  ,k  direction. 

Finally,  we  compare  the  result  (3-72)  to  that  obtained  by 
considering  i     element  of  the  crystal  to  be  a  local  longitudinal 
polarization  modulator.   In  this  approximation  we  have 


Eo  ~?EnRi  sin( *&/*.) 


where  7  is  a  factor  related  to  r  ,  $  =  $(-T)  is  the  crystal  surface 
potential,  and  $  is  the  full-wave  retardation  voltage  which  for  r^ 
modulation  is  given  by 


\   =  X/nl'r63  =  2*/nl3r63ko 


Substitution  for  $  yields 


E2  ~  ?E0R1  Sin  K3r63ko*(-T) 


Then,  using  the  Taylor  expansion 

x3 
sin  x  =  x  -  —  + 
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we  have,   for  small  <J>, 


VE0R1  ~  Inl3r63ko*("T) 


Clearly,  the  first-order  term  above  is  almost  identical  to  the  form  of 
(3.72),  except  for  the  absence  of  the  0  dependence.   It  would  be  of 
interest  to  obtain  an  experimental  evaluation  of  the  resolution  as  a 
function  of  Q    .      If  the  predicted  9     dependence  is  not  observed,  then 
the  first-order  approximation  is  not  adequate,  so  that  the  perturbation 
series  must  then  include  higher-order  terms. 
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SUMMARY 


Our  aim  in  the  present  work  has  been  to  develop  a  theoretical 
analysis  of  the  modulation  characteristics  of  an  electro  -  optic  light 
valve.   The  development  has  been  based  on  Maxwell's  equations  and  a  physical 
configuration  which  is  representative  of  that  of  a  typical  light  valve.  We 
have  attempted  to  use  a  minimal  number  of  assumptions  and  to  justify 
approximations  by  order  -  of  -  magnitude  estimates.   The  fundamental  assump- 
tion is  that  the  linear  electro  -  optic  effect  may  be  expressed  by  an 

r  21 

appropriate  optical  dielectric  tensor. 

In  Section  II  we  obtained  the  exact  expression  for  the  static  field 
present  in  the  anisotropic  light  -  valve  crystal,  in  terms  of  the  Fourier 
composition  of  an  arbitrary  surface  charge  distribution.   In  Section  III  we 
used  the  linear  electro  -  optic  tensor  to  express  the  crystal  wave  equation 
as  a  function  of  the  static  field.   A  first  -  order  perturbation  treatment 
was  used  to  obtain  an  approximation  for  the  optical  fields  transmitted 
through  the  crystal.   This  approximation  was  then  used  to  obtain  an  expression 
for  the  resolution  characteristics  of  an  intensity  -  modulation  light  valve. 
The  development  for  the  solutions  of  both  the  static  fields  and  the  optical 
fields  was  based  on  a  generalized  method  of  images. 
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